The anomalous tunneling of Bose-condensate excitations 
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We discuss the tunneling of phonon excitations across a potential barrier separating two con- 
densate bulks. It is shown that the strong barrier proves to be transparent for the excitations at 
low energy e. Moreover, the transmission is reduced with increasing e in contrast to the standard 
dependence. This anomalous behavior is due to an existence of the quasiresonance interaction. The 
origin of this interaction is a result of the formation of the special well determined by the density 
distribution of condensate in the vicinity of a high barrier. 
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The investigation of the tunneling phenomena in the 
Bose-condensed systems becomes one of the most evolv- 
ing trends in the BEC studies of ultracold atomic gases, 
(see e.g., [1] and references therein). One of the interest- 
ing aspects in this field is associated with the problem of 
the tunneling of collective excitations possessing specific 
features in the condensed phase. The condensate moving 
at velocity v smaller than some critical velocity v c trans- 
mits across a potential barrier without scattering and 
reflection. Excitations of a condensate at rest, in partic- 
ular, sound excitations should experience the tunneling 
transmission together with the reflection. It is signifi- 
cant that the motion of excitations takes place at the 
background of the inhomogeneous built-up distribution 
of condensate density g (x) . The inevitable reduction of 
g (x) near the barrier results in the appearance of a po- 
tential well. For the sufficiently high barrier V ^S> fi 
where \i is the chemical potential, the tunneling of ex- 
citations under these conditions, as we will see below, 
demonstrates an anomalous character. On the one hand, 
the barrier proves to be transparent for phonon excita- 
tions within a limited range of low energies. On the other 
hand, the coefficient of tunneling transmission T reduces 
as the energy of excitations e grows. This is in contrast 
to the typical growth of T with e. Such anomalous be- 
havior to a great extent is associated with the appear- 
ance of a quasi-bound state in the continuous spectrum. 
The formation of the quasi-bound state is connected with 
the specific behavior of function g (x) in the vicinity of 
the barrier. Provided the size L of the condensate is 
much larger than the correlation length £ and the bar- 
rier is high, the dependence g (x) in essential degree has 
an universal character. As a result, anomalous picture 
of tunneling acquires a relatively universal character as 
well. In order to make the exposition ultimate clear, 
we consider one-dimensional symmetrical configuration, 
namely, two identical condensates separated by a rectan- 
gular barrier with height Vo ^> \x. The general equation 
for the Heisenberg field operator of particles ^ (x,t) has 
the form 



n,-i*---nr**. (i) 



Hereafter it is implicitly assumed that the transverse size 
of the system is sufficiently large. Accordingly a collision 
of atoms has a three-dimensional character and, as usual, 
U = 4irh 2 a/m where a is the scattering length. We con- 
sider only the case a > 0. 

We start from the consideration of the ground state, 
assuming temperature T = 0. In this case the operator 

can be replaced with the macroscopic condensate wave 
function W (x) . The equation for 'f (x) in the dimension- 
less form can be written as 



9* 



Id 2 * 



1 )V + V (x)* 



|*| 2 §. 



(2) 



Here we introduced the following notations 

x = x/£, i=nt/h, * = ¥(l7//i)*, V = V//i. (3) 

Further we omit bar for all variables in Eq.(||). 

Let us find a stationary solution (%) of Eq. (j2|) , as- 
suming that (x) — 1 for |x| ^> 1. Outside barrier 
Eq.(0) has the known static solution 



(x) = tanh (\x\ + xq) 



(4) 



Under the barrier we can find the first integral of Eq. (j|) , 



taking into account the condition (0) = 



(5) 



Here b 2 = * 2 (0) , x 2 = 2 (V - 1) + b 2 . The general 
solution takes the form 



*o (x) 



cn 



(6) 



where cnu is the Jacobi elliptic function (see, e.g. Q). 

In the case of strong barrier when xq ^> 1 and x$d > 1 
(d = 2a is the barrier width), one can neglect the nonlin- 
earity related to the interparticle interaction. Equating 
solutions (@),(||) and their derivatives at the boundary of 
the barrier |x| = a, we derive 



(x) = (kq sinh^oo) 1 coshxx, *S>o (a) 



< 1. 



(7) 
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Let us now consider the tunneling of excitations prop- 
agating at the background of the distribution of the con- 
densate density. These excitations can be found as oscil- 
lations of the classical field of a condensate (see, e.g.J|). 
For this purpose, let us represent function ^ in Eq.(||) as 
"F = + ^' an d linearize the equation in Hf' 



i—=h9' + g(x) *' 



(8) 



h = -^ + V(x) + 2g(x)-l, g(x) = H> 2 (x). (9) 

For excitations with energy e, we will seek the solution 
of a set of equations (g) for and 'J 7 '* in the form || 

*' = u(x)e- %et -v* (x)e Z£t . (10) 

Then we have 

(h — e^j u — gv = 0, — gu + (h + ej v = 0. 

Introducing notations 

S = u + v, G — u — v (11) 

and combining these equations, we arrive at fourth-order 
equation 



(h + gj(K- g) 



-q)S = S 2 S. 



(12) 



If the solution of S (x) is known, one can readily find 
function G (x) 



eG=[h-g)S. 



(13) 



Far from the barrier at |x| 3> 1, we have g (x) — > 1. The 
solutions of Eq. ( fL2| ) in this uniform region are sought as 
S ~ exp (ikx) . As a result, we find four roots 



fci, 2 =±fc, k = V2 (VTTe 1 - 1 
^3,4 = Tiq, q = V? (Vl + £ 2 + I 



)1/2 

1/2 



(14) 



The first two roots corresponds to the Bogoliubov spec- 
trum (in the dimensional variables k — ► fc£ and e — » e//i 
respectively). The solutions corresponding to roots £3,4 
do not contribute into the asymptote of the general solu- 
tion at I a; I 3> 1 but prove to be significant in the regions 
close to the barrier. For e *C 1, we have 



k = e, q = 2, 



(15) 



neglecting the terms of the order of e 2 . 

Outside barrier the general solution can be represented 
as a linear superposition of basic solutions Si (x) of equa- 
tion (O) corresponding to the roots given in (114). Let 



excitations fall from the left hand site. Then, omitting 
the divergent components, we have 



S (x) = ASi (x) + BS 2 (x) + CS 3 (x), x< -a 
S (x) = DSi (x) + FS 5 (x) , x > a. 



(16) 



Here the functions S n (x) for fixed e have an universal 
form 

S n (x) = e iknX ( tanh (-x + x ) + % ) , n = 1, 2, 3 



£4,5 (ar) = #2,3 • 

(17) 

Function G (x) can be found from ( fL3| ) with conserving 
the same set of the coefficients to be determined as in 
©■ 

Consider now the underbarrier region. Since g (x) < 
x^ 2 <C 1, we first find the solution of Eq.(JlJ) omitting 
g (x). In this case Eq.([l2|) can be rewritten as 



1£_ 

2dx 2 



V-l] S = e 2 S. 



(18) 



Inserting the solution as S ~ exp (— xx) into the equa- 
tion, we again find four roots 



xi, 2 = ±x+, x 3 ,4 = ±x_, x± = x oy /l ± 2e/x: 2 . 

(19) 

Accordingly, the general solution for the underbarrier re- 
gion reads 

Sb = Ke"+ m + Le-"+ x + Me"- X + Ne~"- X . 



Substituting this expression into ([L3|) and involving ([L9|) , 
we arrive at 

G B = - {Ke*+ X + Le~^ x ) + Me >c - X + Ne' 3< - X . 

Now we can employ that e/xq -C 1 and, supposing 
ea/xo <§; 1, one can replace x± with xq. As a result, 
we have 



S B = (M + if) e^ * + {N + L) e 



(20) 



G B = (Af - K) e* ax + (N-L) e-"° x . 
The boundary conditions at x = ±0 yield eighth equa- 
tions for the determination of all coefficients in ( |l6|) and 

D 

5 (±o) = S B (±0) (a) , G (±a) = G s (±a) (/3) , 



rfS- 

f/.T 



±0 



dx 



(7) 



dG 



dG f 



±0 



(^) 



±0 



(21) 

In order to find the solution in the explicit form, we 
restrict ourselves by considering the tunneling and re- 
flection of collective excitations in the sound region cor- 
responding to the roots (|l5|). In addition, we simplify a 
set of equations omitting quadratic terms like k 2 , fc£ , 



3 



£ 2 (C = <C 1) and involving the smallness of co- 
efficients C and F of the order of (A;, C) with respect 
to the other coefficients in Eq. (|l6|) . This simplification 
yields the equalities M w K and N w L resulting from 
Eq.(pT}5). Together with the boundary conditions for G 
(PH<3) this gives 



2 ' 2 2 



(22) 



Here 



F = Fe- 2a , 
D = De ika , 



B = Be 1 



A = Ae 



— iko 



The remaining set of four equations can be solved 
straightforwardly 



K = 



-k 2 + (k 2 + 4C 2 ) e- A *° a ), D = §4:ik(e- 2 *° a , 



\(k(k + 2iC)e- 3>coa , 



Z = k 2 -(k 



2i( f e - iiiaa . 



(k 2 e~ 



(23) 

As a result, returning to dimensional variables, we find 
for the tunneling transmission of the excitations 



T 



m{k/x ) 2 e - 4 *°a 



(24) 

Even for a wide barrier the transmission coefficient 
demonstrates practically the full transparency at 



1 2 



— 1 a 



For k ^> A:* and still the transmission coefficient 

falls with increasing k 



T 



(25) 



This result is in a drastic contrast with the typical en- 
hancement of T with the growth of k. 

Let us consider a more accurate solution for the un- 
derbarrier region, involving the condensate density g (x) 
in Eqs.(p^|),([l^). The explicit form for g (x) can be ob- 
tained from the results @-(@) for ^ (x) . The numerical 
calculation should be performed in order to find the so- 
lution of Eqs.(12),(|l3j) for the underbarrier region. The 
general solution implies again the use of the boundary 
conditions (|2"T|). The results of calculations for the de- 
pendence T on k are presented in Fig.l. In these figures 
the function T (fc) , obtained with ignoring g (x) in the 
underbarrier region, is shown by the dashed lines. The 
exact solution shows again the existence of the full trans- 
parency. The only difference is a shift of the maximum 
of T to k — 0. With increasing k, in particular, in the 
region where the asymptotic expression ( p5] ) is correct, 
the both curves coincide. We intentionally choose quite 





Figure 1: The coefficient of tunneling transmission T, as a 
function of fc£ for two sets of parameters: (a) xo£ = 8, xod = 
4; (b) >toi = 3, xod = 8. Dashed line corresponds to T(k£) 
obtained with neglecting g(x) in the underbarrier region. 



different sets of the parameters in Fig.l(a,b) in order to 
demonstrate rather universal character of the T (k) be- 
havior. In fact, we see only the change of the scale of 
/c-region where the transparency is high. 

The anomalous behavior of the tunneling transmission 
of excitations has an interesting origin. In the vicin- 
ity of the barrier the condensate density decreases as 
tanh 2 [(|ar| + Xo) /£], producing potential wells for the ex- 
citations on the both sides of a barrier. In the case of 
a strong barrier the ratio (|a;| + xo) /£ <C 1. In these 
conditions a virtual resonance level in the continuous 
spectrum close to e = appears. This results in the 
drastic change of the tunneling picture. As k increases, 
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the quasi-resonance scattering decreases entailing the un- 
usual reduction of T(k). At fe£ ~ 1 the role of virtual 
level reduces and the normal growth of T with k should 
be restored. We can see it from the inserts in Fig.l(a,b) 
where the dependence T (k) is plotted just for this inter- 
val of k. Note that in all cases the anomalous tunneling 
can be revealed considering the reflection coefficient 



R = 



l-T. 



(26) 



Treating the tunneling of phonon excitations, it is in- 
teresting to trace the formation of energy flux Q. In 
the general case this flux can be found employing the lo- 
cal form of the energy conservation law. Provided the 
Hamiltonian is represented as 



H 



E (r) d 3 r, 



h 2 . - U - - - - 

E (r) = — V*+V* + (V - u) * + — 
2m 2 



then 



Using Eqs.@ and (Eh, we find 



Q - -—Re ( 9 * + 9 ^ 
m \ dt dx 



(27) 
(28) 

(29) 



Assuming as before that the condensate is at rest, we con- 
sider the energy transfer by excitations. Treating them 
as oscillations of the classic field of a condensate, we have 
from (p^) for the flux averaged over the time 



<Q)=n ^— (Q), (Q) 



3/2 



111 



1/2 



Re 



dt dx 



(30) 



Here no is the total particle density, Q is the expression 
for the energy flux defined in terms of dimensionless units 
(0) . Let us consider again the incident flux of excitations 
with energy e. Inserting the function \P' from (|Io| ) and 
using notations (|H|), we obtain 



where R is the reflection coefficient (|26J). According to 
( |j~3|) the contribution into ( ]3l] ) from the function G equals 



Im 



G *^D = |A|2 (1 " R) 5 1 1 " tanh2 {x ~ Xo) ^ 



A sum of both expressions is independent of x, demon- 
strating the conservation of the energy flux of excitations. 

Let us find the relation between \A\ 2 and the physical 
parameters. Comparing the density of noncondensate 
particles n' = ^no \A\ 2 in the incident phonon flux with 
the known Bogoliubov result (see e.g. Q) for the relation 
between n' and phonon density n(k) concentrated in the 
fc-mode n' — (n/hcok), we find 



\A\ 2 w 2( A t/Sc fc)(n (k) /n ). 



(32) 



Here c = y/n/m is the speed of sound defined for the 
density at \x\ ^> a. 

Returning to the dimensional units, we find for the 
energy flux (p9|) 



(Q) = n (fc) cos (1 - R) . 



(33) 



In order to obtain the flux under the potential barrier 
in the explicit form, we neglect again the distribution 
of the condensate density in this region. Taking into 
account the solution ( po| ) and the conditions M w K 
and N w L, we find from © (Q) = AexlmAKL*). 
The coefficients K and L are determined in (|23) with 
|A| 2 according to (J 
acquires the value 



As a result, the energy flux p 



(Q) = n (fe) c eT, 



(34) 



where T is the transmission coefficient (p4|). At s > a 
the flux has the same value (|34|), Expressions ( |3^ ) and 
(0) demonstrate the constancy of the energy flux in the 
system (see Eq.(p6|)). Note that this is not the case for 
a particle flux. One can show that this flux inevitably 
connected with the transfer of phonon excitations, is not 
conserved and depends on x. 

This work is supported by INTAS (project 2001- 
2344) , by Netherlands Organization for Science Research 
(NWO), and by the Russian Foundation for Basic Re- 
search. 



(Q> = |lr J) [S*^- + G" d(; 



dx 



dx 



(31) 



Outside the barrier we use the general solution for S 
in the form ( |i"6j ) and ([j"?]). Taking into account that 
function S3 is real and S2 (k, x) — S* (fc, x) = Si (—k, x) , 
we find for the region x < — a 



Im fs*^)=\A\ 2 



k 

R) - (l + tanh 2 (x - x )) 
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